Mobility in a one-dimensional disorder potential 
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In this article the one-dimensional, overdamped motion of a classical particle is considered, which 
is coupled to a thermal bath and is drifting in a quenched disorder potential. The mobility of the 
particle is examined as a function of temperature and driving force acting on the particle. A frame- 
work is presented, which reveals the dependence of mobility on spatial correlations of the disorder 
potential. Mobility is then calculated explicitly for new models of disorder, in particular with spatial 
correlations. It exhibits interesting dynamical phenomena. Most markedly, the temperature depen- 
dence of mobility may deviate qualitatively from Arrhenius formula and a localization transition 
from zero to finite mobility may occur at finite temperature. Examples show a suppression of this 
transition by disorder correlations. 

PACS: 05.40, 05.60, 71.55J 



I. INTRODUCTION 



In recent years there has been wide interest in trans- 
port properties of disordered media: e.g. diffusion on a 
polymer in an external field, random resistor networks, 
domain wall dynamics of magnets in a random field, and 
pinning of vortices in type-II superconductors. Some re- 
views have already been devoted to this subject 

This article focuses on the mobility of a particle moving 
at finite temperature in a one-dimensional disorder po- 
tential. Its purpose is twofold: (i) The functional depen- 
dence of mobility on temperature and on an additional 
external driving force is for the first time worked out ex- 
plicitly in terms of stochastic properties of the disorder 
potential in the continuous space. Previous publications 
have mainly focused on the situation, where the particle 
moves on a lattice. The dynamics was specified in terms 
of hopping rates between neighbored places. Our point 
of view will be more adequate and also physically more 
transparent in situations, where the disorder potential is 
well characterized and hopping rates, if to be used, had 
to be calculated from this potential first. In other works, 
where space was treated as continuous, only special disor- 
der types have been considered ||[| . (ii) This functional 
dependence is explicitly calculated and discussed for new 
models, in particular for spatially correlated distributions 
of the disorder force. Thereby we obtain generalizations 
of the Sinai model JjJ with spatially uncorrelated forces, 
which has attracted particular attention in the past (see 
e.g. H and references therein). 

Our approach, which leads to closed analytic expres- 
sions, is limited to the evaluation of the mean velocity. 



Thus interesting transport phenomena, like an anoma- 
lous scaling behavior of the (mean squared) displace- 
ment as a function of time, which characterize dynamical 
phases Q and were found for the Sinai model [^|^J^,^,|| , 
are beyond the scope of the present treatment. 

In the following, we first derive and discuss the general 
expression for mobility (Sec. |n|). Its asymptotic behav- 
ior for large or small temperatures or driving forces is 
then analyzed (Sec. III). In Sec. [Tv| mobility is calculated 
over the complete parameter range for some models. Fi- 
nally, the phenomena of thermal activation encountered 
thereby are summarized (Sec. fv|) . 



II. BASIC DESCRIPTION 

Our problem is defined by the one-dimensional 
Langevin-equation for a single particle with coordinate 
x in the presence of a disorder potential U(x) and an 
external force F, 



x(t) =F-U'(x(t))+r)(t). 



(1) 



The thermal random force rj(t) is assumed to be Gaus- 
sian distributed with moments 



(#)) = o, 
(r,(t)v(t')) = 2TS(t 



t>). 



(2a) 
(2b) 



Angular brackets represent thermal average in a heat 
bath of temperature T. 

We are interested in the velocity-force-characteristics 
(VFC) for a given disorder potential, where the average 
velocity is defined by 
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v T (F) := lim - (x(t) -x(0)). 

t— >oo t 



(3) 



A bar denotes the translation-average for the given real- 
ization of disorder, i.e. one replaces U(x) by U(x — xo) 
and averages over all xo- In the last formula, this means 
an average over all initial positions x(0). If desired, an 
additional average over an ensemble of disorder realiza- 
tions should be taken as last operation. 

For convenience, we suppose for the moment periodic- 
ity of the potential, U(x + L) = U (x), and formulate the 
result in a way, which does not depend on the periodic- 
ity L. We argue then, that the result is correct for any 
unbiased disorder potential. 

The dynamics of the model can be reformulated in 
terms of the Fokker-Planck equation for the probability 
density P(t,x) and current density J(t,x): 

d t P(t,x) = - d x J(t,x), (4a) 
J(t, x) = [F — U'{x)] P(t, x) — T d x P{t, x). (4b) 

From the stationary solution, which has a homoge- 
neous current distribution, one derives the mean velocity 
v = JL/ J Q L dxP(x) in a standard way (see e.g. fLofl ) 

v T (F) = TL (l-e- LF ' T )/ 

/{ ( L dxe-^' T ( L dx' e ^')/ T + 

y Jo Jx 

+e- LF ' T £ dx e-^/T £ dx , e 0(x')/T | (5) 

with the effective potential cj>(x) :— U(x) — Fx. This ex- 
pression simplifies by replacing the periodicity L by NL 
and taking the limit N — > oo. The result can be repre- 
sented in two ways. A first one, frequently used in the 
literature, will be indicated in the following paragraph. 
The second one, which seemingly has been disregarded 
up to now and which can be evaluated more easily, is the 
basis of the remainder of this article. 
The first version has the compact form 



v T (F) 



-(x) 



where 



t(x) := T~ 



±00 



dx i e [4>(x')-<j>(x)]/T 



(6a) 



(6b) 



with the integral running towards ±oo for F ^ 0. 
This representation of the VFC evaluates the so-called 
sojourn-time density t(x): dt = r{x)dx is proportional 
to the conditional probability of finding the particle in 
the interval (x, x + dx) at times t > 0, provided it started 
from position x at t = 0. Therein all later passages of the 
particle parallel and opposite to the direction of the driv- 
ing force are included. Several works have been devoted 



to the study of sojourn-time distributions, mainly in spa- 
tially discrete hopping models (see e.g. and references 
therein). 

In contrast to Eq. (6), which emphasizes the "dynam- 
ical" aspect of the problem by the statistical analysis of 
the sojourn-time, we prefer a second "static" point of 
view, formulated in terms of correlations of the disorder 
potential. This formulation enables us also to interpret 
our results directly in terms of activation processes in 
the energy landscape of the disorder potential. It is re- 
lated to Eqs. (6) by a mere change of the order of two 
integrations: the mobility 



Ht{F) := v T (F)/F 



can be calculated directly as 



with the generating function 



gr(y) :=exp{[U(x + y)-U(x)]/T}. 



(J) 



(8a) 



(8b) 



Viewed as a function of temperature, gr(y) is the gen- 
erating function of the potential energy difference corre- 
lations at distance y. Eq. (|a|) shows, that mobility is 
essentially the Laplace transform of the generating func- 
tion. In an experimental situation, where the force- and 
temperature dependence of mobility are known, the gen- 
erating function may thus be determined by an inverse 
Laplace transformation. 

The above derivation was based on the fact, that the 
current density is spatially constant in the stationary 
state, which is a particularity of one-dimensional prob- 
lems. Unfortunately this approach does not allow an 
evaluation of the diffusion constant, which requires joint 
probability distributions at different times, and relax- 
ation phenomena, which are not stationary. Both topics 
have been treated for the Sinai model ■ 

Formulae (6) and (8) have been derived for periodic 
potentials. Since periodicity shows up only implicitly as 
a property of t(x) and griy) but no longer explicitly 
as parameter, we postulate their validity also for non- 
periodic disorder potentials. Depending on the nature 
of disorder it may happen, that the stationary current is 
zero for some range of forces. In this case the integral in 
Eq. ( |8a| ) will diverge, leading to vanishing mobility. 

In contrast to the original equation of motion, the ex- 
pressions for the mobility are not invariant under the 
transformation U(x) — > U(x) + const x x and F — > 
F + const. They require an unbiased disorder potential, 
i.e. 



U'(x) = 0. 



(9) 
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This condition, which is obvious for periodic potentials, 
has to be imposed on the non-periodic case as well. How- 
ever, situations with a biased disorder potential can be 
treated, too. Then in our expressions U has to be taken 
as the original potential after subtraction of its bias and 
F as the original external force plus the mean force of 
the original potential. 

For the spatially discrete version of the model with 
uncorrelated hopping rates, it was shown | fi"H that the 
results, obtained from a periodic potential and taking 
the limit of infinite periodicity in the end fs|Jl2|| , are un- 
changed, if one allows for non-periodicity from the very 
beginning. In addition, for zero temperature but arbi- 
trary disorder, one can integrate the equation of motion 
Eq. (0) after separation of variables and finds directly the 
zero temperature limit of Eq. ( |Sa| ) without use of period- 
icity, leading to Eq. (O) below. 

Now we address the question, whether an additional 
average over different realizations of the disorder poten- 
tial may affect the VFC. In principle, this average should 
be taken as last average in Eqs. ( pa] ) or (|a]). Again, for 
the spatially discrete version of the model with uncorre- 
lated hopping rates, it was shown jyj, that velocity is a 
self-averaging quantity. This means, that its value, when 
calculated for a given realization of the disorder poten- 
tial, coincides with probability one with the value after 
an additional averaging over all realizations of the disor- 
der potential. It seems natural to assume this property 
for any disorder potential with short-ranged correlations, 
since the particle samples during its drift the potential 
on infinite length scales, where a given realization is ex- 
pected to be representative for the whole ensemble of 
realizations. Therefore we may consider the spatial aver- 
ages as averages in the ensemble of disorder realizations. 
If, for a contrary example, the disorder extends only over 
a finite region, the ensemble-average does modify the re- 
sult and is to be performed in addition. This situation 
has been studied for the Sinai model without bias . 

Our main formula (8) can be illustrated physically in 
the following way: Consider jumps over a finite distance 
y. The larger the fluctuations of U(x + y) — U(x) when 
x varies, the larger will be gr(y) > 1- The inequality 
holds for all unbiased potentials. Therefore we might call 
Ay(y) := T^ 1 In gx(y) > "energetic roughness" on dis- 
tances y. A large roughness signifies a pronounced relief- 
structure of the potential to be overcome by thermal 
activation and reduces mobility. An "enthalpic rough- 
ness" At(F) can be introduced by exp[AT(F)/T] := 
(T/F) f ±o ° dyexp[-yF/T + A T (y)/T], reminding of the 
relation between free energy and enthalpy in equilibrium 
thermodynamics, since F and y are thermodynamically 
conjugated variables. This enthalpic roughness turns out 
to act as an effective activation energy, since it deter- 
mines mobility by ^t{F) — exp[— A T (F)/T}. 

Eq. (8) shows also, that in the evaluation of the dis- 



order potential an energy scale is set by temperature, 
whereas the spatial structure of the potential is relevant 
only up to the length scale T/F. 



III. LIMITING CASES 

Before we attempt to evaluate the VFC for particular 
models, defined by a probability distribution of the dis- 
order potential, we analyze different limits of the general 
expression for mobility. 



A. Low temperatures 

At strictly zero temperature, we obtain from an inte- 
gration of the equation of motion (|l|) over a finite time 
interval (fcj, i/): 



x f 



t f - U 



= F 



x f 



Jxi 1- 



dx 



(10) 



U'(x)/F 



with Xij = x(tij). If the disorder force is everywhere 
weaker than the external force, U'(x)/F < 1, the particle 
cannot be trapped in the disorder potential. In the limit 
tf — ti — > oo, Eq. ( |l0|) then confirms the zero-temperature 
limit of Eq. (8): 



MF) = {[l-U>(x)/F]-i}~ 



(11) 



In the opposite case, with U'(x)/F > 1 somewhere, the 
particle will be localized, i.e. have vanishing mobility. 
The threshold-forces F c for the onset of drift towards 
x = ±oo are clearly given by the maximum/minimum 
slope of U. In the localized region, Eq. ( |To| ) is invalid in 
the strict sense, since the particle does no longer sample 
the whole potential. However, if it is used naively, local- 
ization formally shows up by a divergence of the integral 
in the denominator. 

At small, but finite temperatures, one might expect 
to find an Arrhenius-like thermally activated behavior. 
This is certainly true for periodic potentials, where one 
easily derives (for F ^ F^ 1 ) 



max ii 



r _ L y -^miH e -[0(Zmax)-0(Zmi„)]/T 



2tt 



(12) 



with A max ,min : — (^max,min) j where £ m in / max de 

notes the position of a minimum/maximum of <f> such, 
that (a; max — x m - m )/F > and the energy difference 
4>( x max) — ^(^min) is maximal. This expression is valid 
only for temperatures much smaller than this activation 
energy and L\F\/T ^> 1, such that activation over the 
maximum occurs only in the direction of F. In this case, 
the mean velocity is just proportional to Kramers tran- 
sition rate (see e.g. Q) for thermal activation out of 
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minimum £ m ; n over x max . The resulting velocity is finite 
for all finite external forces. 

In the case of non-periodic disorder potentials, where 
arbitrarily large energy barriers (or curvatures at ex- 
trema) occur with finite probability density and thus no 
highest energy barrier exists, deviations from Arrhenius- 
like temperature dependence may occur. 



B. High temperatures 

For high temperatures one may expand the generating 
function into 



9t(v) = E "1 T ~" [U(x + y)-U(x)Y 

A ^ Til 



(13) 



An additional temperature dependence of mobility comes 
in through the fact, that the generating function is eval- 
uated at distances y = £,T/F, which increase with tem- 
perature. If qualitatively [U(x + y) — U(x)] n ~ \y\ nCy for 
\y\ — > oo with some roughness exponent £ < 1, we will 
have anyway 



Ht(F) = 1-0 



/ T 2(f-l)^ 



(14) 



On the other hand, if the potential is so rough that £ > 1, 
this series diverges and the particle can be expected to 
be localized at all temperatures. 



C. Small driving forces 

In the limit of small F, the function <?t(?/) is evalu- 
ated in (^a|) mainly for large arguments and the integral 
therein acts as additional spatial average. Then one has 



e U(x)/T e -U(x)/T 



(15) 



provided that the averages therein exist. Therefore, one 
has finite mobility < ji < 1 for small driving forces, i.e. 
Ohmic behavior. The range of this Ohmic regime is given 
by the condition, that T/F has to be small compared to 
typical length scales (periodicity or correlation length) 
of the disorder. Remarkably, mobility is independent of 
the direction of the driving force, even for asymmetric 
potentials. If, on the other hand, one of the averages is 
infinite, mobility vanishes for small forces. In general, 
any non-Ohmic behavior requires the divergence of one 
of the averages in Eq. ( |l5| ) . 

It is interesting to note, that a disorder potential may 
give rise to a temperature dependence of the Ohmic mo- 
bility, which is qualitatively different from Arrhenius-like 
behavior. In order to illustrate this, let us assume a po- 
tential energy density 



P(I7i) := S(U! - U(x)) 
with an asymptotics 

P{U) w cexp[-\U/E\ a } 



(16) 



(17) 



for U — > ±oo with some energy scale E and exponent 
cr > 1. For small temperatures we may evaluate the inte- 
grals in Eq. ( [T5| ) with a saddle-point approximation and 
obtain 



tPr(F) 



cr(cr-l) / E 



2ttc 2 E 2 \<tT 



x exp 



-2(<r-l) 



oT) 



x 



(18) 



Only in the limit a — > oo, i.e. a narrow distribution 
P(U), one obtains an Arrhenius-like temperature depen- 
dence of the velocity. Potentials of a finite width, like pe- 
riodic or quasiperiodic ones, belong to this class. For gen- 
eral a > 1, mobility will vanish with temperature faster 
than Arrhenius-like, since the exponent a /(a — 1) in the 
exponential function will be greater than one. However, 
for a > 1 and T > mobility will always be finite. 

In the marginal case a = 1 we find a transition from 
zero mobility at T < E to finite mobility at T > E. In 
the case a < 1 the particle will always have zero mobility. 

These results show the naively expected tendency: the 
broader the distribution P(U), the larger are the involved 
activation energies and the smaller is the mobility. Per- 
haps less expected is the fact, that spatial correlations 
of the potential for small y will affect the VFC only for 
larger driving forces. 



D. Large driving forces 

In the limit of large F, the function gr{y) is evaluated 
in Eq. ( |3a| ) only for small arguments. We now expand 

9T ^T/F) = 1 + (F±/F) i + \(F±/F 2 ) i 2 + --- (19) 

for C ^ 0, where we introduced in general temperature- 
dependent characteristic forces of order n by 



T-) 9 t(v) 



(20) 



From the representation ( pa[ ) for the mean velocity, we 
obtain 



MF) 



1 



l + F^/F + FZ/F* 



(21) 



If the nature of the disorder potential is "smooth" , i.e. 
allows to change the order of differentiation an averaging, 
one derives from the definition (8b) of g, that 
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U'{x) = 0, 



FJ = TU"(x) + U' 2 (x) >0. 



(22a) 
(22b) 



If in addition U"(x) — 0, as for periodic potentials, 
the leading correction of the asymptotics is temperature- 
independent. 

Otherwise, if the potential is not smooth, the identities 
(22) need not hold, as we will see in the Gaussian model 
below, where may be different from zero. Then dis- 
order affects a finite shift of velocity with respect to the 
disorder-free case even for large velocities, 



v T (F -> ±oo) = F - F± + 



(23) 



for F 2 > F±. 

In the case of a discontinuous potential, one may find 
fl , T(0 ± ) 7^ 1 inhibiting the expansion (|l9|), see e.g. the 
random-potential case of the Gaussian model below. 



IV. EXAMPLES 

We now explicitly calculate the VFC for certain disor- 
der potentials to demonstrate the richness of phenomena 
which can be found in this class of diffusion problems. 

A. Sinusoidal model 

For the sake of completeness, let us discuss in the 
present context the periodic model U{x) — E sin(2Trx/L), 
which has been studied already some time ago by Ambc- 
gaokar and Halperin jl5| in order to determine the ther- 
mal noise contribution to the dc Josephson effect. Af- 
terwards, we will discuss a non-periodic generalization of 
this model. 

The generating function can easily be calculated as 



g T {y)=I a [2{E/T) S m{-Ky/L)]. 



(24) 



where Iq is a modified Bessel function. Due to the 
smoothness and symmetry of the potential, griy) is even 
and analytic. A straightforward evaluation of Eq. j8a| ) 
leads to 



tPr{F) 



sinh Trf/$ 



d£ x 



x cosh[£7r//tf] Jo[(2/i?)cos(£tt/2)] 



(25) 



with reduced temperature i) := T/E and reduced force 
/ := F/Fq where Fq := 2irE/L. This expression can be 
treated further analytically only in limiting cases. The 
result of its numerical integration is shown in Fig. ^. 
At T = one finds from Eq. (0) 



for 
for 



I/I<1 
I/I>1 



(26) 



with a rapid onset of motion near the the critical force 
F c :=rnax x [|tf'(z)|] = F . 

At low temperatures and within the force-region |/| < 
1, where the particle is localized at zero temperature, we 
find from an saddle-point evaluation of Eq. (|25|) 



v T (F) = 2F smh(7r//i?) f 2 e~ EA( - F ^ T (27a) 

E A (F) := 2 J B v /l-/ 2 + 2F/arcsin/, (27b) 



provided $ <C \J\ — f' 2 and ■& <C / arcsin /. The temper- 
ature range of its validity shrinks when the external force 
vanishes or approaches the critical force. This asymp- 
totic evaluation yields the mean velocity just as the dif- 
ference of Kramers rates for an activation parallel and 
antiparallel to the driving force across energy barriers 
EA(F) I f!rfE. Since the energy barriers assume only cer- 
tain values, temperature dependence is Arrhenius-like. 

At high temperatures, we can expand the denominator 
of Eq. ( p5| ) in E/T. This produces only even terms and 
we find to second order 



fPr(F) 



1 



2{E/Tf 



A+{2fE/Tf 



At small forces follows from Eq. (|25| 
(i T (F) « ME/T)]- 2 



(28) 



(29) 



in agreement with a direct calculation from Eq. (|l^) with 



,±u(x)/t = l (E/T). 



(30) 



The characteristic forces of lowest order are i 7 ^ = 
and F^ = Fq /2, since U is smooth. In Fig. [I] one ob- 
serves for finite temperatures, that the VFC approaches 
in the limit F — > oo first the curve T = before it as- 
sumes the limit /i = 1. This reflects the fact, that the 
leading correction to /i = 1 is temperature-independent. 

Now we address the question, how these results will 
change when the periodicity of the potential is aban- 
doned. Let us construct such a potential by stretching 
the maxima U = E from points to plateaus such, that 
the density of minima U = —E is reduced to p < 1/L. 
The generating function griy) behaves for small y as 

g T (y) » (1 - P L) + P LI Q [2{E/T) sin(^/L)] , (31) 

since this distance falls with probability (1 — pL) onto a 
plateau and with probability pL into a trap region with 
U < E, where Eq. (^4|) holds. This expression is correct 
in 0(y 2 ), thus the leading correction for large driving 
forces is given by 



F% = 2pL{nE/Lf 



(32) 
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and vanishes with density p. However, the critical force 
F c = 2irE/L at T = 0, being the maximum slope of U, 
does not depend on p. On the other hand, if there is 
no long-range translational order between the traps, we 
deduce using Eq. (J30|) 



for which one obtains immediately 



e ±u(x)/T „ (1 _ pL ) e ±E/T + pL j Q ( E / T y 

and the asymptotics 

p T {F -> 0) « [2(1 - pL)pLsmh(E/T)I Q (E/T) - 
+ (l- P L) 2 +p 2 L 2 I%(E/T)}- 1 . 

This expression displays a low-temperature behavior 

-l 



(33) 



(34) 



Ht{F -> 0) 



(1 - P L)pL^2irE/T + p 2 L< 
x27T(E/T)e- 2E / T . 



(35) 



Its exponential factor is Arrhenius-like with activation 
energy 2E. The prefactor can not be calculated from 
Kramers expression (12), since the curvature of the max- 
imum is infinite. Remarkably, mobility at zero force and 
its leading correction at large forces depend only on the 
mean density of traps and not on more information about 
their distribution. This is true for any shape of the traps, 
at small forces since exp[±U(x)/T) depends only on the 
shape and density of traps, at large forces since the lead- 
ing term of grin) — 1 for small y is simply proportional 
to the density. 



B. Gaussian model 

In general, the probability distribution V of the func- 
tions U may be generated by an "hamiltonian" H., 



V[U] cx cxp{-H[U}}. 



(36) 



The hamiltonian can have direct physical significance. 
Imagine that the particle diffuses on a substrate. In an 
external field, the shape of this substrate will determine 
the energy U of the particle. Assume that the shape 
of the substrate was frozen after a quench from an ini- 
tial temperature, where the substrate performed shape- 
fluctuations according to a reduced hamiltonian TL (sub- 
strate energy divided by temperature). Then this hamil- 
tonian determines the distribution of potentials U ac- 
cording to Eq. ([36]). 

In this section we consider the Gaussian case, where 
the hamiltonian is bilinear in the potential. We assume 
the "energy" of a realization of the potential to be de- 
termined by a inversion-symmetric stiffness e in Fourier- 
space according to 



H[U} = 



1 



dk 
2^ 



U(k)e(k)U(-k), 



(37) 



9t{v) 



1 f dk 1 — cos(fcy) 



(38) 



The stiffness should satisfy lim^o k 3 /e(k) = = 
lim^oo k/e(k), otherwise the particle is always localized. 
We restrict us to e(k) = mo + mi\k\ +iri2k 2 , thereby gen- 
eralizing Sinai's random-force model with a term rn^fc 2 
only. For simplification we select the cases with m,\ = 
or mo = 0. 

We start with mi — 0, where the absolute fluctuations 
of U{x) are confined. One has 

g T (y) = exp { (l - e"l» ^ yo ) VoFt/t) (39) 

with the length- and force-scales 



2/o := V m 2/mo, 
F T := l/2m 2 T. 



(40a) 
(40b) 



Since the potential now typically is rough on short length 
scales, grin) has a non-analytic distance dependence at 
y = 0. In addition, the first characteristic force 



1% = ±F T 



(41) 



is finite and depends on temperature, as well as the char- 
acteristic forces of higher order. 

Laplace integration over g yields (see Fig. ||a) 



Pt(F) 



1 



(yoF T /T) 



Vo\F\/T p -y F T /T 



yo\F\/T 7(y \F\/T,y F T /T) 



(42) 



where 7 denotes the incomplete gamma function. This 
model exhibits the asymptotic behavior 

Ht(F -> ±00) = 1 - F T /\F\, (43a) 
Ht{F -» 0) = expj-yo^T/T}. (43b) 

It provides an explicit example of a non- Arrhenius-like 
temperature dependence of the mobility, since yoFx/T — 
1/(2^/tooTO2 T 2 ). The fact, that temperature occurs 
squared, is due to the Gaussian nature of disorder in 
agreement with the discussion leading to Eq. ([l8]). Note, 
that the limit F — > always contains the constant m2, 
which enters the distribution 



[ U 2 

P([/)cxexp^- — 



dk 



27re(fc) 
exp{- % /m m 2 U 2 }. 



(44a) 
(44b) 



It would be therefore misleading to reason, that this 
limit should depend only on m2, which governs e(k) for 
small fc. Rather, for this limit all length scales matter, 
as already visible in Eq. (BS) . 
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The special case m± = mi = describes an uncorre- 
lated potential and Eq. ( |l5| ) applies for all forces. As the 
amplitude of fluctuations diverges in the limit m-i — > 0, 
the generating function diverges since lim/^oo k/e(k) ^ 
and the particle is localized at all temperatures and driv- 
ing forces. 

Consider now the case mo = 0, where 



gr(y) =exp{(2 2/1 F T AT)[C , 
-ci |y/yi|cos|y/yi| 



(e c \y/yi\) 



(2 yi F T /nT) 



-si |y/yi|sin|y/yi|]} (45) 
(46) 



with a length-scale y\ :— vnilm\, the force-scale Ft ■= 
\/2miT as before, Eulers constant C = 0.577 . . ., and the 
cosine and sine integrals ci and si. Fig. |2|b shows mobility 
obtained by numerical integration of Eq. (p5|). The ap- 
proximation in Eq. ([46| ) refers to |y/yi| 3> 1. The power- 
law behavior of g for large y translates under Laplace 
transformation to a power-law behavior for small F: 



T[{2 yi F T /^T) 
1 



y/4 Vl F T /T 



[ire 



l - c \F\/2F T f VlFT ^ T) (48) 



r denotes the gamma function, the approximation is 
valid for F <C Ft, the limit refers to T — > 0. Since the 
exponent 2yiFT/nT — l/nmiT 2 is positive, the mobil- 
ity always vanishes with F. The exponent is independent 
of the parameter TO2. Mobility behaves at low temper- 
atures like J' 1 / 7rm i T . l n this case, we observe an even 
stranger deviation from Arrhenius-activation. Since V[U] 
is invariant under U(x) — > U (x) + const, the potential is 
unbound and Eq. ( |l8| ) is violated. 

The case mo = 7712 = again leads to complete local- 
ization, since lim^oo k/e(k) ^ 0. 

The case mo = m\ = is, as already indicated above, 
the Sinai random- force model Q. To show the equiva- 
lence of our approach with the discrete approaches ||[ll| , 
we reproduce the VFC of this case. Eq. (J30|) reduces to 



(49) 



(50) 



3 T (j/)=exp{|y|F T /T} 
for all y and results in (compare Fig. |^c) 

for \F\ < F T 

1 - F T /\F\ for |F| > F T . 

Therefore Ft is the threshold force for the onset of mo- 
tion. The temperature dependence of the model is com- 
pletely contained in this threshold force. In particular, 
one finds at constant F a localization transition at tem- 
perature T c = l/2m,2-F. 

All cases with m2 ^ discussed above have a common 
asymptotics for large driving forces, which is determined 
by 77i2 alone. The reason therefore is, that this asymp- 
totics is governed only by the short-scale fluctuations of 
the potential, determined by leading term of the stiffness 
for large k. 



C. Poissonian model 



Finally, we consider a model for diffusion on a stepped 
crystalline surface, as illustrated in Fig. ||a. This model 
shows asymmetric transport properties and also exhibits 
a localization transition. 

The profile of the surface is described by the height 
function h(x), which increases in units of a > for in- 
creasing x. The structure is characterized by the prob- 
ability distribution S(l) for the step length I. We use 
a Poisson-distribution S(l) = exp(—l/lo)/lo with mean 
value Iq. The particle is supposed to be in constant 
force- field with components f Xj h in direction parallel and 
perpendicular to the a;-direction, which might be com- 
ponents of a gravitational or electric field. They give 
rise to the total potential 4>{x) — f%h(x) — f x x. In or- 
der to use our previous formulae, we identify U (x) — 
f h h(x) - (a/l )f h x and F = f x — (a/l )f h by subtracting 
the mean bias. 

From the Poissonian distributed number of steps in an 
interval of length y one obtains easily 



g T (y) = exp[j/F±(T)/T] (51) 



for ^ ^ with asymmetric forces of first order 



F±(T) = -{a/l Q )f h T (1 - e ±a ^ T )T/l % 0, (52) 



which depend on temperature. From Eq. 
(see Fig. |b) 



lh we obtain 





1-F±(T)/F 



for 
for 



F^T) < F < F+(T), 
F % F±(T). 

(53) 



Since at T = the localization region covers —(a/lo)fh < 
F < oo (corresponding to < f x < oo) and shrinks to 
zero for T — > oo, we find for forces in this region a local- 
ization transition at a temperature T c , which is implicitly 
determined by F — F^(T), as depicted in Fig. ||c. 

Apart from the asymmetry, the Poissonian model is 
very similar to the random-force Gaussian model. In- 
deed, the force is distributed independently at different 
sites, this time, however, with an asymmetric and singu- 
lar distribution. 
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V. CONCLUSION 



An analysis of particle mobility in one-dimensional dis- 
order was presented. We developed a framework based on 
the generating function of spatial correlations of the dis- 
order potential. For arbitrary temperature and driving 
force, these correlations are relevant up to a length T/F. 
At large driving forces, mobility depends only on local 
properties of disorder, whereas at small driving forces 
global aspects matter. 

For some models, which generalize previously studied 
structures of disorder, mobility was evaluated over the 
complete range of temperature and force. Thermally ac- 
tivated motion led to a rich phenomenology 

The temperature-dependence of mobility can deviate 
drastically from Arrhenius formula. This is character- 
istic for systems with a broad distribution of activation 
energies. Such deviations have been obtained already by 
an ad hoc averaging of escape times over a spectrum of 
activation energies (see e.g. |l6|]). This procedure is un- 
satisfying from a principal point of view, since transport 
properties depend also on the spatial location and not 
only on the height of energy barriers, which is in higher 
dimensions even more important than in one dimension. 

The models with spatially uncorrelated force distribu- 
tion, the Poissonian model as well as the random-force 
Gaussian model, exhibit a localization transition. The 
similar phenomenology of these models is due to the va- 
lidity of the Central Limit Theorem, which assures a 
Gaussian distribution of U(x + y) — U(x). 

It is therefore of particular interest to examine the sta- 
bility of the localization phases with respect to the intro- 
duction of correlations in the force-distribution. For the 
Gaussian model, we achieved this by two different, con- 
tinuous deformations of the random-force model. Switch- 
ing on the couplings mo or mi, we found the localized 
phase to be unstable, i.e. mobility was finite for all posi- 
tive temperatures and driving forces. 

This result complements previous studies: In a dynam- 
ics, where the particle may hop on a lattice only in one 
direction and which is thus incompatible with a Langevin 
equation of motion, dynamical phases were found to be 
stable with respect to short-ranged correlations in the 
hopping rates, but unstable with respect to long-ranged 
correlations p7| . Introducing force-correlations into the 
Sinai model, the scaling behavior of diffusion in the ab- 
sence of a driving fore has been found to be modified, 
but remained anomalous [ ^8|Jl9| ]. 
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FIG. 2. Mobility of the Gaussian model as a function of 
the driving force for different temperatures T = 4, 2, 1, 
0.5, 0.25, 0.125 from top to bottom in each diagram: (a) for 
mo = m2 = 1 according to Eq. (^), (b) for mi = mi = 1 af- 
ter numerical Laplace integration of Eq. (^), and (c) for the 
random force model as limiting case mo = mi = 0,m2 = 1. 
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FIG. 3. Poissonian model: (a) sketch of geometry, (b) mo- 
bility according to Eq. ( |53| ) as a function of the driving force 
at temperatures T = 4, 2, 1, 0.5 from top to bottom and 
(c) localization temperature as a function of driving force for 
afh =lo = l. 
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